Domain Wall Dynamics in Brane World and Non-singular Cosmological Models 
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We study brane cosmology as 4D (4-dimensional) domain wall dynamics in 5D bulk spacetime. 
For a generic 5D bulk with 3D maximal symmetry, we derive the equation of motion of a domain 
wall and find that it depends on mass function of the bulk spacetime and the energy-momentum 
conservation in a domain wall is affected by a lapse function in the bulk. Especially, for a bulk 
spacetime with non-trivial lapse function, energy of matter field on the domain wall goes out or 
comes in from the bulk spacetime. Applying our result to the case with SU(2) gauge bulk field, 
we obtain a singularity-free universe in brane world scenario, that is, not only a big bang initial 
singularity of the brane is avoided but also a singularity in a 5D bulk does not exist. 



I. INTRODUCTION 

In the standard model of cosmology, we usually ex- 
pect an initial singularity in the beginning of the uni- 
verse. It is the so-called big-bang singularity. In fact, 
Penrose-Hawking singularity theorems [1] prove that our 
universe must start with a singularity. If a singularity ex- 
ists, physical quantities diverge, and then classical gen- 
eral relativity will break down there. We need the so- 
called "quantum gravity" , which is respected to resolve 
a singularity problem, however it has not been found. 
Although quantum gravity is one of the most important 
subjects in physics, we are now not only looking for quan- 
tum gravity but also trying to construct a unified theory 
of all fundamental interactions. 

As such a unified theory, a string/M theory is pro- 
posed, which is expected to include quantum gravita- 
tional theory [2,3]. It suggests that we are living in a 
3-brane in a higher dimensional spacetime [4,5]. In such 
a brane world picture, particles in the standard model 
are expected to be confined to 3-dimensional (3D) brane, 
whereas gravitons propagate in the entire bulk spacetime. 
As a result, the extra dimensions can be rather large 
(<0.1 mm) [6,7]. An alternative brane model is also pro- 
posed [8] . It has been found that 4D gravity is recovered 
on the brane in low energy limit. Based on such a new 
world picture, a lot of brane world cosmological mod- 
els have been extensively studied and found new aspects 
in cosmology [9-12]. Since the spacetime is inhomoge- 
neous even if the 3-brane is maximally symmetric just 
as a Friedmann- Robertson- Walker spacetime, the analy- 
sis will be complicated. First, the Friedmann equation 
was obtained by solving 5D inhomogeneous spacetime, 
and two new effects including dark radiation term were 
found [9]. Then new method to obtain covariantly the 
effective 4D gravitational equation was found [10]. An- 
other way to study such a brane cosmology is a domain 
wall dynamics. The motion of a domain wall in a given 
spacetime provides the Friedmann equation of a brane 
world [11]. It was found that the mass of Schwarzschild- 
anti de Sitter (AdS) background black hole spacetime 



corresponds to dark radiation through its tidal force on 
a domain wall. 

Although a brane scenario provides new aspects in cos- 
mology, it does not resolve a singularity problem. There 
is a big bang singularity in a simple brane cosmological 
model based one the Randall-Sundrum model. Except 
for the case with no dark radiation, a black hole singu- 
larity also exists in the bulk although it is hidden beyond 
the horizon. 

As for a singularity problem, if a string/M theory re- 
ally gives a unification of all fundamental interactions 
and particles, it will be free from a singularity. Then we 
expect that the effective theory will show such a singu- 
larity avoidance. The effective string/M theory predicts 
that a bulk space may contain many fields, such as a 
dilaton scalar field or gauge fields. In particular, if our 
5D bulk spacetime is effectively obtained by dimensional 
reduction of the 11D Hofava-Witten model [3], a dilaton 
scalar field and U(l) gauge field appear [13]. Also in an 
effective action approach of a string theory, the Gauss- 
Bonnet term as the leading order quantum correction 
may appear to gravity, in particular, in the case of het- 
erotic string [14]. Also using brane structure, new mech- 
anisms of spontaneous symmetry breaking of gauge field 
have been proposed [15]. For example, the present stan- 
dard model, that is SU(3)xSU(2)xU(l) gauge theory, is 
obtained on the brane in higher dimensional SU(5) gauge 
theory. It suggests that higher dimensional bulk contains 
non-Abelian gauge field in the brane world scenario. In 
this viewpoint, a bulk may not be simply Schwarzschild- 
AdS spacetime, but it would be a various type of ge- 
ometries by the effects of bulk fields. Therefore, in this 
paper we investigate how geometry of a bulk spacetime 
contributes to dynamics of domain wall. 

In the case with U(l) gauge field in a bulk, it is shown 
that an initial singularity is avoided [16-18]. However, 
since the bounce of a domain wall occurs always inside of 
the inner horizon of Reissner-Nordstrom-AdS black hole, 
such a bounce may not be realistic because the inner 
horizon may be unstable and then it will evolve into a 
singularity [19]. In this paper, we study dynamics of a 
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domain wall in a spacetime containing non-Abelian gauge 
field in the bulk, and try to resolve a singularity problem. 
This is because there is a globally regular spacetime, so- 
called particle-like solution, as well as non-trivial black 
holes in the 5D Einstein-SU(2) Yang-Mills system with 
spherical symmetry [20]. We adopt SU(2) Yang-Mills 
field as non-Abelian gauge field. 

In the following section, we first derive the equation 
of motion of 3D domain wall in the 5D bulk space- 
time. In Sec. Ill, after reviewing a brane cosmology in 
the pure Randall-Sundrum case and the case containing 
U(l) gauge field, we discuss the model with SU(2) bulk 
gauge field. We find a singularity-free brane cosmological 
model, in which both a bulk singularity and a big bang 
initial singularity are avoided. Summary and discussion 
follow in Sec. IV. In appendix, we discuss black holes and 
a particle-like solution of Einstein- Yang-Mills system. 



II. DYNAMICS OF DOMAIN WALL 

The evolution of the universe in brane cosmology is de- 
scribed by the dynamics of 3D domain wall (our universe) 
moving in a 5D bulk spacetime. It strongly depends on 
the spacetime structure of bulk structure. Hence we in- 
vestigate the domain wall dynamics in a generic bulk 
spacetime with 3D maximal symmetry. We assume that 
the 3D domain wall T> is homogeneous and isotropic, 
and a bulk spacetime is static. Just as in the Randall- 
Sundrum model [8], we assume that a 5D cosmological 
constant A5 is negative. We introduce a typical scale 
length Z by Z 2 = — 6Ag . Assuming the existence of 3D 
maximal symmetry, the general form of 5D metric with 
a negative cosmological constant is written as 
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k takes the value of 0, 1, or —1, corresponding to flat, 
positive, or negative curvature of 3D maximally symmet- 
ric space, which is the same geometry of domain wall T>. 
The bulk geometry is represented by two functions, that 
is a lapse function 8{r) and a mass function p(r). In this 
coordinate, a domain wall is located at r — a(r). We 
assume the Z 2 -symmetry between the two sides of the 
domain wall. 

The equation of motion of domain wall is given by 
Israel's junction condition [21] at the domain wall, 



curvature of T>, where is unit normal vector of T>. 
T^u is the energy-momentum tensor of 4D matter fields 
on the domain wall. In Gaussian normal coordinate, we 
rewrite the metric (2.1) as 



ds 2 = dr] 2 + ry ab dx a dx b 



(2.4) 



near the domain wall T>, where '-fab the induced metric in 
this coordinate ststems. The condition (2.3) becomes 
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where 



K ab = -n^dyTlab- 



(2.5) 
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Assuming Z 2 -symmetry (K+ b = —K~ b ), we glue two 
patches of the bulk spacetime which boundaries are given 
by a domain wall at r = a(r), where r is a proper time 
on the domain wall [22]. In the given bulk metric (2.1), 
the 5-velocity of the domain wall = dx^/dr and the 
unit normal vector n M are defined by 
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where a dot denotes the derivative with respect to r. 
The minus sign appears because the unit normal vector 
is taken to be the inward direction. Then we evaluate the 
extrinsic curvature from Eq. (2.6). Their spatial compo- 
nents are given by 
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Generalizing the calculation in [23], we find the tt com- 
ponent as 



K+ T {= -K~ T ) 
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Now we take the energy-momentum tensor of a domain 
wall in the conventional form as 



Tab = PU a U b + P{jab + U a U h ). 



(2.11) 



This perfect-fluid form includes normal matter fluid de- 
scribed by p m and P m as well as a domain-wall tension 
a, i.e. p — p m + a and P = P rn — a. Using Eq. (2.9), 
spatial part of the junction condition (2.5) is reduced to 
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where h^ v — — n^n^ is the induced metric of the 
hypersurface V and — h°hffl is the extrinsic 
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where H = d/a is the Hubble expansion rate, and the 
rest gives 
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Eq. (2.12) with Eq. (2.2) gives the equation of the do- 
main wall, 
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where 
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Eq. (2.14) corresponds to the Friedmann equation of our 
brane universe. When p(a) is a constant, i.e. the bulk is 
the Schwarzschild-AdS spacetime, we recover the conven- 
tional brane universe model, which contains a quadratic 
term of energy density and dark radiation [11]. 

Here we should stress one point about the energy- 
momentum conservation law of matter fluid on the brane. 
Taking a derivative of Eq. (2.12) and comparing it with 
Eq. (2.13), we find that 



p + 3H(P + p) = ^pd. 

da 



(2.17) 



This means that the energy-momentum of the matter on 
the brane is not conserved unless the lapse function 6(a) 
is constant. We discuss two possibilities in order. 
(1) Suppose that the tension vanishes and the equation 
of state is given by P m = (7 — l)p m . Then we find the 
equation for p m as 



p m + 3-t Hp m = Sp mi 
which is integrated as 
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(2) In the case that the tension does not vanish but is 
constant, we find 



■ 3-yHp m = S(p rn + a), 
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the integration gives 
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The second term proportional to a tension appears fur- 
ther in addition to non-trivial contribution of 6(a). If 
the 5D bulk is asymptotically AdS spacetime, then 8(a) 
converges to such that energy conservation law is recov- 
ered as the universe expands. But it is just a sufficient 
condition. In fact, it turns out in next section that it 
is realized even in quasi-asymptotically AdS bulk space- 
time. In the next section, taking a concrete model, we 
shall discuss this term in detail. 



III. DYNAMICS OF BRANE UNIVERSE 

In order to analyze the motion of a domain wall, we 
rewrite the Friedmann equation (2.14) in the form of a 
particle motion with zero energy in a given potential as 



\a 2 + U(a) 



0, 



(3.1) 



where the potential U (a) is given by 



U(a) = - 



p(a) A 
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where p m {a) is given by Eq. (2.21). The expansion law is 
modified by a bulk matter flow into (or from) the brane, 
in addition to quadratic energy density and 'dark radia- 
tion'. In what follows, we discuss the Friedmann equation 
for several brane models in order. 



where ag is an initial scale factor, Sq = S(ag) and po is an 
integration constant, corresponding to the initial energy 
density. Since p m oc a~ 37 in standard cosmology, the en- 
ergy in a proper volume will decrease when the universe 
expands if the lapse function S(a) is a decreasing function 
with respect to a. However, if 5 converges to constant as 
a — > 00, then the energy is asymptotically conserved as 
pm — * Goa~ 37 . This case includes the equation of state 
for vacuum energy (P m = —p m ), which does not keep 
constant in the evolution of the universe. Then if we 
tune the parameters in order to find zero cosmological 
constant initially (A 4 + 87rG 4 po = 0), the present small 
vacuum energy might be explained by the difference be- 
tween initial and final values of vacuum energy, although 
the lapse function must be an increasing function. 



A. Vacuum Bulk Spacetime 

If bulk spacetime is vacuum just as Randall-Sundrum 
model, the background is described by Schwarzschild- 
AdS spacetime [24], that is p(r) = M and S(r) = 0. 
Then the potential (3.2) is 



M A 4 2 
k — —a 
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This case was discussed in detail in [11]. Unless M is 
negative, we cannot avoid a big bang singularity. If we 
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assume that there exists an event horizon (otherwise, un- 
known effect from naked singularity will come into our 
brane universe), M is positive for k — 1 or 0. However, 
if k = — 1, M could be negative M > — Z 2 /4. In such a 
case, there is a chance to avoid a big bang singularity. In 
fact, if there is no matter fluid, a big bang singularity is 
avoided. However, when matter fluid appears, if 7 > 2/3, 
the potential diverges to —00 as r — > 0. Then, a big bang 
singularity may not be avoided. We show a schematic 
diagram in Fig. 1. Furthermore, there is a black hole 
singularity in 5D bulk spacetime (i.e. at r = 0). 



B. Bulk U(l) Gauge Field Model 

In this case, the 5D bulk spacetime is a Reissner- 
Nordstrom-AdS spacetime, whose metric is fi(r) = M — 
Q 2 /r 2 and S(r) = 0. It has a timelike singularity at 
r = 0, and two horizons (an event and inner horizons). 
Since the lapse function vanishes, the energy-momentum 
on the brane is conserved. 

The potential (3.2) becomes [16-18] 



U(a) = 
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New term corresponding to U(l) charge Q behaves just 
as stiff matter energy density with the opposite sign. It 
works as a repulsive force in the domain wall dynamics. 
If there is no matter fluid, a solution turns out to be 
a catenary type bouncing universe model when A4 van- 
ishes. The bounce of domain wall occurs near the inner 
horizon. Furthermore, an oscillating universe is realized 
when k = 1. In this model, an initial big bang singularity 
is avoided. 

However, we can show that any bounce occurs inside 
the inner horizon, which is expected to be unstable. In 
fact, we find that /(a) > U(a) when A4 vanishes. The 
bounce occurs at the inner point of U(a) = 0, while the 
horizons exist at the points of f(a) = 0. Then we find 
that the bounce point appears inside of the inner horizon 
[16-18]. We show one typical diagram in Fig. 2. The in- 
ner horizon may not appear in the realistic situation and 
it will evolve into a singularity because of its instability 
[19]. Hence the bouncing universe is not realized in the 
case with U(l) charge. 

We also have to stress one more important point. Once 
we include matter fluid in our model, if 7 > 1 including 
radiation fluid, the potential becomes unbounded from 
below. We naturally find a big bang singularity at r = 0. 
Even for the dust fluid 7 = 1, in order to obtain a finite 
potential, there is a constraint for Q 2 as 



Q 2 > 
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C. Bulk SU(2) Gauge Field Model 

Now we apply our domain wall dynamics to a model 
with SU(2) bulk gauge field. The bulk geometry is ob- 
tained as a solution of 5D Einstein- Yang-Mills system, 
whose action is 



S 



1 



(R - 2A 5 ) 
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Tr(F 



ym) 
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where F YM is field strength of 5D SU(2) Yang-Mills field 
and g is its gauge coupling constant. 

There are two differences from the models in the pre- 
vious subsections. First, there are new types of solutions 
[20]; a globally regular spacetime, so-called a particle- 
like solution, and non-trivial black holes, so-called colored 
black hole solutions, with spherically symmetry (fc = 1). 
Those are numerical solutions and have non-trivial lapse 
function. We can show that such solutions are stable 
against spherical linear perturbations if A5 is negative. 
They are 5D version of the solutions formed in the con- 
ventional 4D Einstein equations by [25-34]. In the case 
of k = 0, we find that there is no particle- like solution but 
colored black hole solutions exist. However, in the case 
of k = —1, there is neither a particle-like solution nor a 
colored black hole solution. There is only analytic black 
hole solution for k = — 1 (see the detail in appendix). 
We should stress that a particle-like solution has no sin- 
gularity, thus such a bulk spacetime is singularity-free. 
Secondly, the obtained solutions arc not exactly asymp- 
totically AdS, following the definition by [35,36]. In this 
bulk spacetime, some conserved quantities such as the 
"mass" diverge. However, any curvature is finite and the 
metric approaches AdS as r — > 00. So this geometry can 
be called quasi-asymptotically AdS [20,37]. In this quasi- 
asymptotically AdS spacetime, a lapse function converges 
to zero as r — > 00, so the energy conservation law of mat- 
ter fluid on the domain wall is recovered at large radius. 
We shall discuss several bulk models in order. 



1. A particle-like spacetime 

First, we consider a brane dynamics in the case that 
bulk spacetime described by a particle-like solution (k = 
1). The bulk spacetime is parametrized by b — w"(r = 
0), where w(r) is the SU(2) gauge potential. The regular 
solution is found for 6 m ; n < b < 0. 6 m i n depends on 
A5 and is obtained numerically (see Table I). The bulk 
spacetime is stable against spherical perturbations in the 
range (6 m i n <) b s < b < 0, where b s is the lower bound 
for existence of a stable solution. Assuming A4 = and 
H ^> p m , we find that a domain wall can be located in 
this bulk spacetime if 6 m ; n < b < 6 cr i t . 6 C rit is maximal 
value such that U (a) across zero, so no timelike domain 
wall exists for 6 cr ;t < b. To find a stable singularity-free 
universe model (b s < 6 cr ; t ), the 5D cosmological constant 
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A5 is also constrained as l/X < 1.199 (Table I), where 
A = (G^/g 2 ) 1 ^ 2 which is a typical scale length in the 
Einstein- Yang-Mills system. Since 



G 5 = 



GaI 



(3.7) 



which is obtained by Eqs. (2.15) and (2.16), assuming 
A4 = 0, the restriction for a gauge coupling constant is 
given by 



g 2 < 0.719 



(3.8) 



We depict a typical potential and bulk geometry (metric) 
in Fig. 3. Thus, we can realize a singularity-free brane 
universe in the range b s < b < 6 cr ;t, i.e. both big bang sin- 
gularity and singularity in the 5D bulk spacetime can be 
avoided. From the potential in Fig. 3, we easily find that 
a domain wall dynamics leads to an oscillating universe 
model. For the case A 4 < 0, oscillating solutions are 
obtained similarly, but the range of b becomes smaller. 
For the case A4 > 0, catenary type universe solutions are 
obtained for all particle-like bulk solutions. 



2. A black hole spacetime 

Next, we consider the case of a bulk black hole space- 
time. Since we know the behaviours for trivial black holes 
with 6(r) = 0, we discuss only non-trivial ones, i.e. col- 
ored black holes. In this case, a geometry of domain wall 
allows k = 1 or 0. (In the case of k = — 1, a bulk space- 
time can be the analytic black hole solution, so domain 
wall solutions are similar to those in Reissner-Nordstrom- 
AdS case.) Such black hole solutions are classified into 
two types by existence of an inner horizon. If an inner 
horizon exists, we find a bounce solution in domain wall 
dynamics. We show a typical potential in Fig. 4 and bulk 
geometry in Figs. 5,6. In fact, a catenary type universe 
can be realized for k = 0, while an oscillating universe 
can be obtained for k = 1 . Just as the case of Reissner- 
Nordstrom-AdS bulk spacetime, the domain wall always 
enters inside of the event horizon and escapes through the 
inner horizon to another quasi-asymptotically AdS space- 
time. However, an inner horizon may be unstable, then 
a singularity will appear on the Cauchy surface. Such a 
brane universe is not a good candidate for a singularity- 
free cosmological model. 

If a bulk spacetime does not have any inner horizon, a 
domain wall must start from an initial singularity. In this 
case, because of a bizarre behaviour of interior structure 
of colored black hole, U(a) have infinite relative maxi- 
mums at the point on which /z(a) approaches most, so 
domain wall expands repeating acceleration and decel- 
eration. We show a typical potential in Fig. 7 and bulk 
geometry in Figs. 8,9. In the case of k = 1, a domain wall 
expands to a maximum radius a max , where a max satisfies 
U(a max ) = 0. After then, it contracts to a singularity 



repeating acceleration and deceleration. In the case of 
k = 0, a domain wall continues expanding to infinity. 

We summarize the behaviours of a domain wall in Ta- 
ble. II. 



3. Energy-momentum conservation 

Now we examine the energy conservation law of matter 
fluid on the domain wall. From our numerical analysis of 
for a particle-like spacetime, we find that a lapse function 
is almost constant in the stage of decelerating expansion 
(or accelerating contraction). Then the energy in this 
stage is conserved. In the region of accelerating expan- 
sion (or decelerating contraction) , a lapse function is not 
constant, so energy conservation on the domain wall is 
broken. In fact, energy on the domain wall escapes to 
a bulk spacetime as the universe is expanding. For the 
colored black hole bulk, similarly, a lapse function is ap- 
proximately constant outside the event horizon, and the 
energy-momentum is conserved there. Inside the event 
horizon, energy on the domain wall is not conserved and 
escapes to a bulk as expanding. Particularly, much en- 
ergy flows out to bulk at the point which domain wall 
dynamics changes from acceleration to deceleration (for 
the case of colored black hole bulk without inner horizon) 
or at bounce (for one with inner horizon). It shows that 
energy out-flow to the bulk plays an essential role for the 
acceleration of a domain wall. 



4- Asymptotic behaviour 

At last, we examine asymptotic behaviour of a domain 
wall as a — > 00 (or when a is enough large) . The first term 
of Eq. (2.21) describes matter fluid in standard Fried- 
mann equation, which is affected by a lapse function. 
The second term is induced energy through some inter- 
action between a brane tension a and a bulk spacetime. 
Inserting 



5{a y = --(«/ (a)) 2 , 
and the asymptotic behaviour of w as a — > 00, i.e. 



(3.9) 



w (a) -> 



(3.10) 



with a constant K , which was obtained in the appendix 
of [20] , we can evaluate the asymptotic behavior of second 
term as 



2aK$ 
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+ 



C 



3(7 - 2) a 6 a 3 7 
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for 7^2, where C is an integration constant. This first 
term behaves just as a stiff matter, while the second term 
in renormalizcd into the first term in Eq. (2.21) because 
5 — > 0. For 7 = 2, we find 
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-2a K. 



2 ln(a/ai) 



(3.12) 



which behaves similar to a stiff matter fluid but with 
small difference (the existence of In a term). Hence the 
second term of Eq. (2.21) is not so important. 
The 'dark radiation' term behaves 



H(a) 



M n 
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(3.13) 



The first term corresponds to the usual dark radiation, 
while the second term will be more dominant than the 
first term because of In a term, which is still less impor- 
tant than dust matter fluid. 



IV. SUMMARY AND DISCUSSION 

Assuming homogeneity and isotropy for 3D domain 
wall, we investigate the dynamics of domain wall in a 
5D bulk spacetime. The bulk spacetime is described 
by a lapse function S(r) and mass function /x(r). From 
the Israel's junction conditions, we derive two equations 
describing dynamics of domain wall, one is Fricdmann 
equation which is same form to [11] replacing the mass 
of Schwarzschild-AdS black hole by n{r). Another is the 
modified energy conservation law of matter fluid on the 
domain wall. It is worthy of notice that when a lapse 
function S(r) is non-trivial, energy of matter fluid on the 
domain wall is not conserved. 

As for the case with non-trivial lapse function, we 
study 5D non-trivial black hole spacetime called colored 
black hole and globally regular spacetime [20] , which are 
solutions with non-Abelian bulk gauge field. In this case, 
energy of the matter fluid on the domain wall leak out 
to the bulk spacetime by the effect of non-triviality of 
lapse function as universe expands. The energy out-flow 
plays an essential role for the acceleration of a domain 
wall at the accelerating expanding region. On the other 
hand, energy conservation law is asymptotically recov- 
ered at large radius of a. While, the mass function works 
as a repulsive force near a = 0, so a domain wall solution 
with appropriate parameters of a bulk spacetime gives 
an oscillating or a bouncing universe model. It is wor- 
thy of notice that with a globally regular bulk spacetime, 
a singularity-free stable cosmological model in a brane 
world scenario can be realized. Both a big bang singu- 
larity and a bulk singularity are avoided. We summarize 
our result of domain wall dynamics in Table. II. 

In the present analysis, we have studied only SU(2) 
bulk gauge field. In the 4D case, however, because we 
have other gauge symmetries, there are many non-trivial 
solutions with non-Abelian gauge fields [38] . Then we ex- 
pect that the similar solutions with different gauge sym- 
metries exist for 5D case as well. Such a brane picture 
containing non-Abelian gauge fields in the bulk, through 
a dimensional reduction of a unified theory, could re- 
solve the singularity problem. In the effective action of 



a string theory, we also expect that not only standard 
forms of gravity (the Einstcin-Hilbcrt action) and Yang- 
Mills fields exist but also higher-order terms may appear 
in higher dimensions. In higher dimensional Einstein- 
Yang-Mills system with higher order terms of Yang-Mills 
fields, particle-like solutions with finite energy are found 
[39]. It is important to investigate whether in such a 
string theory or its effective theory, a singularity will not 
really appear. 
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APPENDIX A: BULK SOLUTIONS 

In this section, as a non-trivial bulk spacetime, we 
present static solutions of 5D Einstein- Yang-Mills sys- 
tem with a 3D maximal symmetry, i.e. a spherical (k — 
1), plane (k = 0), or hyperbolic (k = —1) symme- 
try. As for spherically symmetric solutions, we discuss 
in detail in [20], where we found a colored black hole 
and a particle-like solution. These correspond to the 
Bartnik-McKinnon solution [25,30] and a colored black 
hole [26,30] in 4 dimensions. Furthermore, an analytic 
solution of non-trivial black hole is found in 5 dimen- 
sions. 

In this appendix, we extend our analysis to the case 
with k 7^ 1. We normalize all physical variables by a 
length scale by A = (G5/.9 2 ) 1 / 2 just as in the previous 
paper. 

We solve the Einstein equations whose action is given 
by Eq. (3.6), assuming the metric form by Eq. (2.1). The 
magnetic part of the vector potential A = A^dx^ of the 
SU(2) gauge field having spherically, plane, or hyperbolic 
symmetry can be obtained by the similar discussion in 
Appendix A of [20], and is described by one function 
w = w(r) as 



A? = 0, A a r =0, A$ = (p,0,w) 
A a n 



{wS k {ip), — — ,0), 



A a 



where 



dip 

( ^(VO s - m Q — w S k (ib) sin 8, cos 9), 

dip 



(Al) 
(A2) 

(A3) 
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simp k = 1 
ip k = 
sinh-0 k = — 1 



(A4) 



and the field strength F = F^dx^ A dx v is written by 



9^A V — d v A^ [Ap, A v 



(A5) 



Note that pure gauge field (F = 0) can be obtained by 
w = ±1 in the case of k = 1 and u> = in the case of 
k = 0, but there is no pure gauge field in the case of 
k = -l. 

With the above ansatz, the Einstein equations and 
Yang-Mills equation are reduced to a set of equations; 



p! = 2r 



f k w' 2 + 



(k-w 2 ) 



2\2 I 



and 



6' = --w' 2 , 
r 



-(rf k e- 5 w'y + \e- 5 w{k - w 2 ) = 0. 



(A6) 
(A7) 

(A8) 



Those with k = 1 correspond to Eqs. (3.1)-(3.3) of [20]. 
In [20] , we discuss a particle- like solution and exterior of 
a colored black hole in the case of k = 1 . Here we analyze 
the other cases (k = and —1) as well as interior of a 
colored black hole with k = 1. 

First, we search for analytic solutions. In the case of 
k = 1, it was found in [20] that there are two analytic so- 
lutions, one is Schwarzschild-AdS spacetime and another 
is non-trivial black hole spacetime. For k = 0, the trivial 
solution corresponding Schwarzschild-AdS spacetime [24] 
is obtained as 



/i = M, (5 = 0, w = 0, 



(A9) 



while there are no non-trivial analytic solutions. In the 
case of k = —1, the absence of pure gauge field yields 
that a magnetic part of SU(2) gauge field cannot provide 
any trivial solution. However there is non-trivial analytic 
solution, as 



p = M + 2\nr, 5 = 0, w = 0. 



Next we search for numerical black hole solutions. Sim- 
ilarly to the case of k = 1, there are two free parameters, 
which are a horizon radius, r^, and a value of potential at 
horizon, % = w(rh)- We find that a black hole solution 
exists in the parameter range < Wh < 1 and any value 
of rh for k = 0. However, for k = —1 we cannot find 
any black hole solution with a magnetic part of SU(2) 
gauge field. We show our numerical results in Figs. 5-11. 
The asymptotic structure at spacelike infinity is similar 
to that in the case of k = 1 , discussed in [20] . The colored 
black hole solutions are classified into two types by exis- 
tence of an inner horizon. In 4D, for the solutions with- 
out any inner horizon the interior of the colored black 



hole shows a very bizarre structure [33,34]. In 5D, those 
black holes also show a similar behaviour. As r decreases 
below horizon, fi decreases and approaches zero, but be- 
fore reaching zero, /U suddenly grows exponentially and 
reaches some very large value. After then, [i falls down 
again to very small value which is closer to zero than be- 
fore. This behavior is repeated until r = (see Fig. 8). 
For a solution with an inner horizon, /i across zero and 
fk begins to increase to very large value after decreasing 
to a certain value (see Fig. 5). 

Finally we have analyzed a particle-like solution. For 
the case of k = and k = — 1 , there is no solutions which 
satisfies regularity at the origin. Then no particle-like 
solution exists. Hence a particle-like solution exists only 
for the case of k = 1. 

Now we discuss stability of the bulk spacetime ob- 
tained above. For asymptotically flat spacetime in 4D, 
Volkov et al. proved that a particle-like solution has n 
unstable modes if the gauge potential w(r) has n nodes 
[28]. It is also true in asymptotically AdS particle- like 
solutions [30]. Since the particle-like solution for asymp- 
totically flat spacetime must have finite node, they are 
unstable. On the other hand there is solutions with no 
node for asymptotically AdS particle-like solutions, so it 
becomes stable. Hence, in 5D case we also expect that 
our bulk spacetime is stable if w(r) has no node. 

To analyze the stability, we perturb the metric and 
gauge potential as 



S(r,t) = S (r)+S 1 (r)e iut , 
w(r, t) = w (r) + wi(r)e luJt . 



The perturbation equations are derived as 

,.2 



Mi 



2r 



- —fJ-1 ^ Wl 



4 

8 1 = 
and 



(All) 
(A12) 
(A13) 



(A14) 
(A15) 
(A16) 



(AlO) _^ (r/oe -^y_ /oe -*v ( i 1 + Sl \ 
r A h \ r Jo J 

+-{rfoe- s °w[)' + i-= ^u;i = -u 2 — 



Wi = -U) -j-lUl, 

Jo 

(A17) 

where fio,5o, wo are given by the unperturbed bulk solu- 
tion and /o = 1 — Mo/?" 2 + r 2 /l 2 . As boundary conditions 
for a particle-like solution, we impose w\ — and w[ = 
at the origin r = 0, because of regularity. The asymptotic 
form of Eq. (A17) are 



r 2 w'{ + 3rw[ = 0, 
which is solved as 



(A18) 
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W!=A+^, (A19) 

where A and B are integration constants. We impose 
A = because of convergence of total energy flux of 
gauge field. We solve them as an eigenvalue problem, 
restricting uo 2 e R. If uj 2 is always positive, the solution 
is stable. We show our result in Table. I. We find that 
if b s < b < 0, uj 2 is always positive, and the solution 
is stable. As we expect, the range of b s < b < corre- 
sponds to those of solutions with no node. Hence we may 
conclude that the solution is always stable if wq has no 
node, while it is unstable if wq has at least one node. We 
also find that for I < 0.779A, there is not any solutions 
with a node, so all solutions are stable. We list up our 
numerical result in Table. I. 
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FIG. 1. The potential U(a) for the case of a vacuum bulk 
spacetime. The solid, dotted, and dashed lines depict U(a) 
with M = 1 for k = 1, 0, — 1, respectively. We also show the 
case with M = — 1 for k = — 1 by a dot-dashed line. 



FIG. 3. The potential 17(a) in a particle-like bulk space- 
time. The solid line depicts the potential U(a). The dotted 
and dashed lines depict a lapse function 5(a) and mass func- 
tion (J,(a), respectively. We set the parameters I = 0.1 and 
b= -7.3. 
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FIG. 2. The potential U(a) for the case with a bulk U(l) 
gauge field. The solid, dotted, and dashed lines depict U(a) 
with M = 1 and Q 2 = 0.1 for k — 1, 0, and —1, respectively. 



TABLE I. The parameter range for which a solution exist. 
The bulk solutions are found in the range 6 m i n < b < 0. They 
are stable if b s < b < 0. A domain wall can be placed in the 
bulk in the range fe m i n < b < fe cr i t . We then find that a domain 
wall can be placed in a stable bulk spacetime if / < 1.199A. 
All bulk solutions are stable for I < 0.779A. 



l/X 




6s 


bait 


10.0 


-0.644 


-0.013 


-0.633 


5.0 


-0.668 


-0.053 


-0.639 


3.0 


-0.718 


-0.148 


-0.661 


2.0 


-0.801 


-0.329 


-0.707 


1.5 


-0.896 


-0.570 


-0.768 


1.3 


-0.958 


-0.735 


-0.810 


1.2 


-0.998 


-0.837 


-0.837 


1.0 


-1.105 


-1.066 


-0.914 


0.8 


-1.178 


-1.270 


-1.035 


0.7 


-1.390 


bmin 


-1.125 


0.6 


-1.552 


bmin 


-1.248 


0.5 


-1.781 


bmin 


-1.423 


0.4 


-2.128 


6min 


-1.691 


0.3 


-2.711 


bmin 


-2.142 


0.2 


-3.883 


6min 


-3.054 


0.1 


-7.411 


6min 


-5.808 
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FIG. 4. The potentials U(a) for a colored black hole space- 
time with an innner horizon. The solid and dotted line depicts 
the potential for k = 1, «Jh = 0.701 and k = 0, tt>h = 0.56. 
We set the parameters as I = 1, r"h = 1. 
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FIG. 6. The lapse functions 5{r) for a colored black hole 
spacetime with an innner horizon. The solid and dotted line 
depicts the potential for k = 1, u% = 0.701 and k = 0, 
TOh = 0.56. We set the parameters as I = 1, rh = 1. 




FIG. 5. The mass functions /x(r) for a colored black hole 
spacetime with an innner horizon. The solid and dotted line 
depicts the potential for k = 1, = 0.701 and k = 0, 
TOh = 0.56. We set the parameters as i = 1, rh = 1. 



FIG. 7. The potentials [/ (a) for a colored black hole space- 
time with no innner horizon. The solid and dotted line depicts 
the potential for k = 1, Wh = 0.4 and k — 0, Wh = 0.7. We 
set the parameters as I = 1, = 1. 
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FIG. 8. The mass functions fi(r) for a colored black hole 
spacetime with no innner horizon. The solid and dotted line 
depicts the potential for k — 1, mi, = 0.4 and k = 0, Wh. = 0.7. 
We set the parameters as I = 1, rh = 1. 



FIG. 10. The gauge potential w(r) for a colored black hole 
spacetime with an innner horizon. The solid and dotted line 
depicts the potential for k = 1, = 0.701 and k = 0, 
w h = 0.56. We set the parameters as I = 1, r h = 1. 
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FIG. 9. The lapse functions S(r) for a colored black hole 
spacetime with no innner horizon. The solid and dotted line 
depicts the potential for k = 1, u% = 0.4 and k = 0, Wh — 0.7. 
We set the parameters as I = 1, Th = 1. 
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FIG. 11. The gauge potential w(r) for a colored black hole 
spacetime with no innner horizon. The solid and dotted line 
depicts the potential for k = 1, Wh = 0.4 and k = 0, wi, = 0.7. 
We set the parameters as I = 1, Th = 1- 
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TABLE II. The behaviour of domain wall for several bulk 
spacetimes. EH and IH denote an event and an inner hori- 
zon, respectively. BB and BC denote the universe which be- 
gins from a big bang singularity, and that ends up with a 
big crunch singularity, respectively. Then, BB^BC means 
evolutionary path of the universe. O denotes an oscillating 
universe, while C means a catenary type universe. 



bulk field 


k 


bulk spacetime 


horizon 


evolution 




1 






BB->BC 


vacuum 


~0~ 


Sch-AdS 


EH 


BB 




-1 






BB 


U(l) or 


1 









SU(2) 


TT~ 


RN-AdS 


EH + IH 


C 


(electric) 


-l 






C 






analytic BH 


EH + IH 







l 


colored BH I 


EH 


BB^BC 






colored BH II 


EH + IH 





SU(2) 




particle-like 


None 





(magnetic) 





colored BH I 


EH 


BB 






colored BH II 


EH + IH 


C 




-l 


analytic BH 


EH + IH 
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